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Abstract 

We use the matrix model to describe the Af = 2 SO(N)/Sp(N) supersymmetric gauge 
theories with massive hypermultiplets in the fundamental representation. By taking the tree 
level superpotential perturbation made of a polynomial of a scalar chiral multiplet, the effective 
action for the eigenvalues of chiral multiplet can be obtained. By varying this action with 
respect to an eigenvalue, a loop equation is obtained. By analyzing this equation, we derive 
the Seiberg-Witten curve within the context of matrix model. 



1 Introduction 



Recently Dijkgraaf and Vafa [1] have made a conjecture, the exact superpotential and gauge 
couplings for a class of Af = 1 gauge theories can be obtained by calculating perturbative 
computations in a matrix model in which the superpotential of the gauge theory is interpreted 
as an ordinary potential. The earlier works [2, 3, 4] motivated this conjecture. Based on this 
observation, there are many works on this direction [5]- [32]. In particular, we restrict to the 
supersymetric SO(N)/ Sp(N) gauge theories. The model with quartic tree level superpoten- 
tial for adjoint chiral field was found [33] and the effective superpotential was computed in 
the context of matrix model and string theory on Calabi-Yau geometry with flux. The per- 
turbative calculation for glueball superpotential was studied in [34]. For arbitrary tree level 
superpotentials, the planar and leading nonplanar contributions were derived by using higher 
genus loop equations and diagrammatics [35]. A field theoretic derivation of the superpotential 
was given in [36] based on the factorization property of Seiberg-Witten curve. An equivalence 
of M = 1 gauge theories deformed from M = 2 by the addition of superpotential terms was 
studied with flavors [37] and without flavors [38] based on the Cachazo-Vafa's idea [2]: The 
low energy information is given by extremization of the effective superpotential and in the field 
theory analysis it is given by characterizing to the factorization locus of Seiberg-Witten curve 
and the equivalence of two description was given in [2]. 

In this paper, we compute the matrix path integral over tree level superpotential obtained 
from M = 2 SQCD by taking arbitrary polynomial of a scalar chiral multiplet as a perturbation. 
The effective theory action can be expressed as a function of an eigenvalue of chiral multiplet. 
The saddle point equation implies an algebraic equation defined on a hyperelliptic Riemann 
surface. The presence of this curve allows us to study the relation of matrix model and the 
gauge theory result of perturbative calculation. By using the basic idea of matrix model, we 
calculate a partition function in terms of a glueball field, a distribution of eigenvalue of chiral 
multiplet, a perturbed superpotential and the mass of quarks. By reading off the two free 
energy contributions from a partition function, one obtains the final effective superpotential in 
terms of homology basis. By varying this effective superpotential with respect to the coefficient 
function appearing in the algebraic curve, one realizes the existence of a meromorphic function 
on Riemann surface with the appropriate structure of zeros and poles. In doing this, the correct 
counting of the number of physical D5-branes in the presence of orientifold planes (05-planes) 
is very important because these values determine the structure of zeros and poles precisely. 
By identifying this function with the resolvent of matrix model, the Seiberg-Witten curves for 
M = 2 SO(N)/Sp(N) gauge theory with Nf hypermultiplets are rederived. For U(N) gauge 
theory with Nf flavors of quarks in the fundamental representation, the derivation of Seiberg- 
Witten curve was found in [39]. It would be interesting to study the results in [40, 41, 42] to 
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deal with the model given this paper. 



2 SO (AT) matrix model 

We will derive the Seiberg-Witten curve for Af = 2 SO(N) gauge theory with Nf fundamental 
hypermultiplets by computing the matrix path integral using the saddle point method [43, 39]. 
Let us consider an Af = 2 supersymmetric SO(N) gauge theory with Nf flavors of quarks 
Q l a (i = 1, 2, • • • , 2Nf, a = 1,2, ■ ■ ■ , N) in the vector (fundamental) representation [44, 45, 46, 
47, 48, 49, 50, 51]. In terms of Af = 1 superfields, Af = 2 vector multiplet consists of a field 
strength chiral multiplet W£ b and a scalar chiral multiplet $ afe both in the adjoint representation 
of the gauge group. The Af = 2 superpotential takes the form 

W tree ($, Q) = V2Ql<S> ab QiJ tJ + V2m tJ QlQi (2.1) 

where Jy is the symplectic metric (® 1 l)®1-N f xN f used to raise and lower SO(N) flavor indices ( 
l NfXNf is the NfXNf identity matrix ) and is a quark mass matrix ( 1 p)(8)diag(mi, • • • , m Nf ). 
Classically, the global symmetries are the flavor symmetry Sp(2Nf) and U(1)r x SU(2) r chiral 
R-symmetry. When Nf < N — 2, the theory is asymptotically free and generates dynamically 
a strong coupling scale A N=2 . The instanton factor is proportional to A. 2 ^r 2 A ~ 2Nf . Then U(1) R 
symmetry is anomalous and broken down to a discrete Z 2 N-iN f -A symmetry by instanton. By 
taking a tree level superpotential perturbation AW made out of the adjoint field in the vector 
multiplet to the Af = 2 superpotential (2.1), the Af = 2 supersymmetry can be broken to Af = 1 
supersymmetry. That is, 

[f] a 

W = W tree (<Z>, Q) + AW, AW = Y, Tr$2fe - 

fc=i 

Then a microscopic Af = 1 SO(N) gauge theory is obtained from Af = 2 SO(N) Yang-Mills 
theory perturbed by AW. 

By using the perturbed superpotential in addition to the tree level one, substituting the 
whole superpotential into the SO(N) matrix model at large iV and replacing the gauge theory 
fields with matrices, we study the various contributions to the free energy. Then the partition 
function can be written as 
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[d$PQ]exp 



vol(SO(N)) 
where for simplicity we change the notation 



-W(^) - V^Qi^QiJij - j2m l3 Q\Qi 
9s 



(2.2) 



W($)= ]T ^Tr$ 2fc . 



k=i 
2 



A superpotential W of order 2( y +1) breaks the gauge symmetry down to a direct product 



of 



+ 1) subgroup. One can write the derivative of W with respect to the field 

N N 

W'(x) = x 2N + E s 2l x^ = n (V - e 2 ) (2.3) 

i=l i=l 

where e^'s are the classical moduli and the symmetric polynomial s 2 k in ef is 



S2k = (-l) k E 



2 2 

e,-, • • • e,- 



»i<— <ifc 

The description for the addition of the mass term for the adjoint scalar only was studied in 
the context of matrix model [52]. The matrix description for pure flavors without any adjoint 
fields to check the Seiberg duality was observed in [53]. Let us study SO(N) matrix model by 
considering even N and odd N case separately because the Jacobian has different form in each 
case and also the spectral curves are different. Now we first analyze SO(2N) matrix model. 
• SO(2N) matrix model 

According to the procedure [3, 54, 39, 35] and by integrating over Q in our case, the 
eigenvalue basis provides 

Z ~ / II W II (>« - X l) ex P ~f E 2^(A a ) - £ log (Xl - mf) 

J o=l a<b " s a=l i=l 

where ±iA a are the eigenvalues of $ = (^q) <S> diag(Ai, • • • , Ajv) and rrii is a quark mass 
(io) ® diag(mx, • • • ,rriN f ). Note the factor 2 in the first term of the exponent. The second 
term comes from the determinant of (<3> + m). The new thing in our problem is the flavor part 
in the last term. For SO(2N) theory without any flavors these terms are absent [35]. After 
exponentiating, the effective action for the eigenvalues is given by 

N i N N f 

S(\) = - E log (A 2 - \l) + - E 2^(A a ) + E log (\l - m\ 

a<b " s a=l i=l 

In this way, the potential W(X) contains a collection of N variables Ai, A 2 , • • • , Aat. Remember 
that due to the antisymmetric property of $ and the trace of it vanishes, only even terms 
in the potential W which is a polynomial of order (2iV + 2) contribute. The saddle point 
equations (classical equations of motion) coming from varying the action with respect to a 
single eigenvalue A a are 

N 9X 1 N f \ 

b^a A a ~ \ 9s i=1 \ ~ m i 

To solve this let us introduce the trace of the resolvent of the matrix $ [43, 3, 54, 39, 35] 

^)4^4£aff- (2 - 5) 

a=l a 
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Then multiplying (2.4) by 2A a / \x 2 -X 2 a ) and summing over an index a, one gets a loop equation 
for uj{x) 

"\x) + lu{x)W'(x) + ^ = 

where the S is defined as S = g s N being fixed in the large N limit and the polynomial f(x) is 
given by 

o=l A a X 

which is a polynomial of order (2iV — 2) with even powers. Therefore the function f(x) deter- 
mines the solution of the matrix integral. Here we take the large N limit and drop the terms 
like uj{x)/x and uj'{x) which will be important when we expand it with respect to 1/N in order 
to derive the Seiberg-Witten differential completely within the framework of the matrix model. 
The spectral curve reduces to 

y 2 = W'(xf-f(x), f(x)=Y:b 2n x 2n 

where we define 



y( X ) = W'(x)+Sw(x). 

This is nothing but a hyperelliptic curve in (x, y) plane. We have to determine the iV unknown 
coefficients b 2n - As in [2], there exists two particular points denoted by P and Q located at the 
two pre images of 00 of x. The force equation becomes 

dS 

MX) = -g. w 
The solution for resolvent is given by [43] 



u(x) = ^W'(xf - f(x) - W\x). 

which is expressed as an iV unknown coefficient function appearing in the polynomial f(x). 
The resolvent has the branch cuts among which the eigenvalues of the matrix are distributed. 
In the large iV limit, the distribution of eigenvalues can be written as 

1 N r 
P(A) = ^£<KA-A a ), j p(\)d\ = l 



a=l 

and the resolvent becomes in this limit 



A — x 



which implies that 



P(A) = ^- + it) - u(X - it)} = ^- [y(X + it) - y(X - it)} . 



The filling fractions are given by 



Si = 



2ni J a 



ydx, 



4iri JAo 



ydx 



where we take the half of the cycle around A due to the orientifold projection. In order to find 
out the functional behavior of f(x) within the matrix model, the saddle point computation of 
the partition function gives rise to up to l/g s term as follows: 



cxp 



+ exp 



-4 / d\p(X)W(X) + ^ / dXdX'p(X)p(X') log(A 2 - A' 5 
- / dXp(X) log(A 2 - m 2 ) 







= exp 


9-^2 J~ 1 




K 9s 



(2.6) 



To get the effective superpotential, one should know both the variation of Ti under a small 
change in Si and T\ that can be read off from (2.6). For the former, we take the following 
change in p(A) 

P(A) - P(X) + 6 -^5(X - a) 

where is an arbitrary point along the i-th cut on the hyperelliptic Riemann surface. From 
the explicit form of JF 2 in (2.6) one considers 



5F 2 = 5S t (W( ei ) - 2S j dXp(X) log(A 2 - e 2 ) 



up to Si independent terms which are not relevant. Then the partial derivative of JF 2 with 
respect to Si becomes 

dJ-'o ~ f p . - .„/■,. ,^( p xdx 



A 2 - x 2 



= -2 [ dxW'(x) -AS f dXp(X) f 
= -2 [ P dx(W'(x) + Stu(x)) = -2 [ P y(x)dx = - ( [ P + [)y(x)dx (2.7) 

Ja Jei \Jei JQ J 

up to an irrelevant constant of integration terms. Moreover the T\ term can be written as 



p xdx 



N f fP 

/ u(x)dx 

Jrrii 



N f 

2S^fdXp(X) 

N f f P 1 N f f P 1 N f f P 

y(x)dx = -J2 y( x )dx + -J2 y(x)dx 

_i Jrrii A Jrrii & -_i J —rrn 



N 



(2.8) 
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up to the Si independent terms. 

Combining the two contributions (2.7) and (2.8) one gets the effective superpotential 



I , P i N f , P i N f r p 

W = --(2N-2) ydx+-^ ydx + -Y, ydx + --- (2.9) 

Z, JQ Z i=1 Jrrii Z i _ l J-m i 

Here there are some remarks in order. In the type IIB string theory, the U (N) gauge theory 
is realized by the worldvolume of N D5-branes wrapped on S 2 and in the dual geometry D5- 
branes are replaced by RR fluxes generating the effective superpotential and the S 2 by S 3 . In 
order to deal with the gauge groups SO(N) and Sp(N), one needs to introduce the orientifold 
plane into the geometry. This will change the contributions of RR fluxes below. The physical 
D5-brane charge of orientifold plane (05-plane) is —1 and the total 2N D5-branes wrapping 
around the origin should be modified by (2N — 2). Moreover the branch cuts in figure 1 in 
[38] are symmetric, due to the Z 2 symmetry, with the one located in the center (We follow the 
notations given in [38]). This implies that the contribution from the compact cycle with a_ k is 
exactly the same as the one with a k . That is, by replacing the D5-branes with the fluxes, one 
gets [38] 

[ h=\N k = I h= l -N^ I h = \{2N -2). 

J Z J d—k Z J clq Z 

By summing over the all a k contour, 



j p h = ( 2 2 N * + 2N o) ~ 2 = ( 2iV ~ 2 ) ~ N f> f h = -(2N -2), £ 



h = -l 

where h is an one form on the Riemann surface. So h should have a pole of order 1 at P 
and Q with residues (2iV — 2 — Nf) and — (2iV — 2) respectively [37]. By using the properties 
C_ fc = — Y!j=i P-j + Co and C k = Z)jLi Pj + Co given in [38], one can divide into two parts: the 
cycle around Co and the one around = 1, 2, • • • , k. Then we have extra two contributions 
in (2.9) denoted by • • •, Yang-Mills coupling term and the term with the cycle (3j. But if we 
take the variation of W , the first contribution will give rise to a trivial cycle and the second one 
gives an element of the period lattice. For U (N) gauge theory, the discussion on this matter 
was considered in [2]. Since Si can be determined by f(x) (and therefore &2n), we have to 
compute the variation of W with respect to b 2n . 

Oy x^ n 

-dx = dx, n — 0, 1, • • • , N — 1 



<%2n 2y 

which are basis for the subspace of holomorphic differentials (one forms) which are odd under 
the Z 2 transformation x — > —x, on the Riemann surface. Note that the full space of holomorphic 
differentials has the dimension (2N — 1) due to the genus (2iV — 1) of the Riemann surface. 
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By changing the bases to the homology basis, the extremum condition of W will give rise 

to 

fQ rP / rrrii r-m,i\ 

(2N-2) Ck-(2N-2N f -2) Cfc - E ( / +/ K fc = 

Jpo Jpo j =1 \Jpo Jpo / 

modulo the period lattice and po is an arbitrary generic point on the Riemann surface. 

There exists a function on the Riemann surface with an (2iV — 2)-th order pole at Q, an 
(2iV — 2Nf — 2)-th order zero at P and simple zeros at x — ±mj for each % = 1, 2, • • • , Nf, 
according to Abel's theorem [2, 39]. The function is simply related to the resolvent divided by 
x 2 (We refer to [37] with flavors and [38] without flavors for the geometric picture) 



z{x) = -*- — = i 2N f <2N - 2 - 

This function has an (2iV — 2)-th order pole at Q and at least fourth order pole at P since f(x) 
is a polynomial of at most (2iV — 2)-th order. Therefore f(x) should contain a factor (x 2 — mf) 
for each % in order for z(x) to have a simple zero at x — ±mj and should be at most (2iVy + 4)-th 
order in order for z(x) to satisfy an (2N — 2Nf — 2)-th order zero at P. Combining these two 
conditions and putting the proportional constant as A 4 ^™ 1 4 , the spectral curve is 



N N. 



f 



2 TT/ 2 2\2 A 4JV-2iVf-4 4 TT / 2 2\ 

y = ll(x -ej -A N=2 x H(x -mj 
i=i i=i 

which is exactly the Seiberg-Witten curve [55, 47, 48, 50, 51] in the field theory analysis where 
we used the relation (2.3). So far we assumed that the number of flavors are not too large. 
Next we consider the case of large flavors. 

For 2N — 2 < 2Nf, the function z(x) is given by 




*W = V -IT- - ~ -±T> 2N-2<2N f <4N-4 



where A(x) is an 2A^-th order polynomial and g(x) is proportional to x^Yl^J^x 2 — mf). The 
presence of this new function comes from the fact that when the 2Nf is greater than (2N — 2), it 
is not enough to have the function f(x) only because the pole structure at P needs to introduce 
a new function. Now we take the proportional constant as A^fl 2 2Ar/ 4 . Then the function z(x) 
vanishes at x — ±mj for each % = 1, 2, • • • , Nf and has an (2iV — 2)-th order pole at Q and 
(2Nf — 2N + 2)-th order pole at P. The expression inside of square root in z(x) should be 
proportional to y(x) 2 /x A 



A (xf A 4iV-2V / -4 2] _ W'{xf 

Jb • -i Jb Jb 



m 
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where f(x) is a polynomial of order at most (2N — 2). The solution for this up to G(A 4 ^ f =2 2Nf 4 ) 
is given by 



N 



.4N-2N f -4 [ 4 
A N=2 



I 

i=l 



AT 



x 4 n(^ 2 -^)-2P(x)n(^ 2 -^) 



i=l 



where P(x) is defined as a polynomial of degree (2Nj — AN + 4) in a; and to*. Therefore the 
spectral curve and the function are given by 



N 



y 2 = m- 2 -e 2 r-m 



i=i 

N 



i=l 



N f 

I< 

i=l 



y A(x) 

z(x) = — — . 

x A x A 

Therefore in both regions of the number of flavors the spectral curve coming from the matrix 
model calculations coincides with precisely the known Seiberg-Witten curve. 
• SO(2N + 1) matrix model 

Let us describe the odd case. Since the presentation looks similar to the one in previous 
discussion, we will present the main difference only. By integrating over Q, the eigenvalue basis 
provides [3, 54, 39, 35] 



r N N „ N i N N f 

Z ~ / J] [d\] II ( X l - X l) II ^ exp -I £ 2W(X a ) - E log (XI - ml) 



N N n N i N 

a=l a<b a=l i> s a=l 

where ±i\ a are the eigenvalues of <3> = ( _ y J ) <g> diag(Ai, • • • , Xn, 0) and m, is a quark mass 
(JJ) ® diag(mi, • • • , rriN f ). Note that there exists an extra factor IT^Li X 2 a in this case [35]. 
Then the effective action for the eigenvalues is given by 

N „ N -i N 

J 1 

a<b a=l V s a=l 

By varying the action with respect to an eigenvalue one gets 



N „ N I N N f 

S(X) = - E log (X 2 a -Xl) - E log X 2 a + — J2 ™(X a ) + E log (Xl - m 2 ) . 

a<b a=l "s a= \ i = i 



f^a _ Xl X a g, 



N f 



i A^ - mf 



0. 



In the large N limit, the second extra term comparing with the previous case SO(2N) matrix 
model does not contribute. Therefore all the arguments from (2.5) to (2.8) in even SO(2N) 
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gauge theory are valid for S0(2N + 1) matrix model. By an appropriate counting the physical 
D5-brane charge of the orientifold which is equal to — |, one gets 

/ h=\N k = [ h= l -N^ [ h = \(2N -\). 

By summing over the all contour, 

j p h = (^ 2 Yl N k + 2N o^ - 1 = (2iV- 1) -N f , jf /i = -(2JV-l), £ h = -l 

where h is an one form on the Riemann surface. 

The extremum condition of W (2.9) where the coefficient (2iV — 2) is replaced with (2iV — 1) 
will give rise to 

rQ fP j!^ / rm,i f-mi\ 

(2N-1) Ck-(2N-2N f -l) Cfc-E / +/ K* = ° 
modulo the period lattice. Then the function is simply related to the resolvent divided by x 



, s V W '( X ) W'ixf fix) w'ix) Ar , 

z(x) = - — = \ V^-^V — , 2N f <2N-l. 2.10 

x x V x 2 x 2 x 

This function has an (2N — l)-th order pole at Q and at least third order pole at P since f(x) 
is a polynomial of at most (2N — 2)-th order. Therefore /(re) should contain a factor (x 2 — mf) 
for each i in order for z(x) to have a simple zero at x = im, and should be at most (2iVy + 2)-th 
order in order for z(x) to satisfy an (2N — 2Nj — l)-th order zero at P. Combining these two 
conditions and putting the proportional constant as A^f^ 2 2iV/ 2 , the spectral curve is 

N N f 
2 TT/ 2 2\2 A 4AT-2iV/-2 2TT/ 2 2\ 

y =[[( x - e J - A w=2 a: H(x - m J 

i=l i=l 

which is known as the Seiberg-Witten curve [56, 47, 48, 50, 51] in the perturbative calculation 
in gauge theory side. 

For 2N — 1 < 2Nf, the function is given by 



, x A(x) 2 g(x) A(x) 

z(x ) = J-LJ--M1J. 2AT - 1< 2iV f < AN - 2 2.11 

V x 2 x 2 x 

where A(x) is an 2iV-th order polynomial and g(x) is proportional to x 2 Y\^l 1 {x 2 — mf). Now 
we take the proportional constant as A^ 2 2JV/ 2 . Then the function z(x) vanishes at x — ±mj 
for each % — 1, 2, • • • , iVj and has an (2iV — l)-th order pole at Q and {2Nf — 2N + l)-th order 
pole at P. The expression inside of square root in z(x) should be proportional to y(x) 2 /x 2 

A ( X Y K AN-2N f -2^ ( 2 2 , _ W'jxf f (x) 

—r~ A ^=2 m J - -2 —^r 

Jb • i Jo Ju 



where f(x) is a polynomial of order at most (2N—2). It is straightforward to see the coincidence 
of the Seiberg-Witten curve. From the properties of z(x) in (2.10) and (2.11), one can easily 
see that 

h(x)dx = — 

z 

is a meromorphic differential with simple poles at P, Q and x = ±ra; with residues (2N — 2Nj — 
2), — (2iV — 2) and 1 respectively for SO(2N) gauge theory, for example. 

3 Sp(N) matrix model 

In this section, we continue to study the matrix model for the symplectic group Sp(N). Let 
us consider an Af = 2 supersymmetric Sp(N) gauge theory with Nf flavors of quarks Q l a (i = 
1, 2, • • • , 2Nf, a = 1, 2, • • • , 2iV) in the fundamental representation. The tree level superpotential 
of the theory is obtained from [44, 45, 46, 47, 49, 57] 

W tree ($, Q) = V2Q l a § a b QlJ bc + V2m tJ QlQiJ ab (3.1) 

where J a b is the symplectic metric ( 1 ~ ) (g) InxN and quark mass matrix (_ ) ® 

diag(mi, • • • , mjvy)- Classically, the global symmetries are the flavor symmetry 0(2Nf) and 
U(1)r x SU(2)r chiral R-symmetry. When Nf < 2N + 2, the theory is asymptotically free 
and generates dynamically a strong coupling scale Aat=2- The instanton factor is proportional 
to A^f^ 2 Nf . Then U(1)r symmetry is anomalous and broken down to a discrete Z 2 N-N f +2 
symmetry by instanton. 

According to the procedure [3, 54, 39, 35] and by integrating over Q in our case, the 
eigenvalue basis provides 



r N 9 N 1 iv 

z - / n w n ( x i - ^i) n ^ exp -± ^ 2w(\ a ) + e ^ - 

J a =l a<b a=l Q " — 1 



N N f 

m 2 



9s a= i j=i 

1 o 



where ±iA a are the eigenvalues of $ = ( _ 1 ) ® diag(Ai, • • • , Ajv) and is a quark mass 
(-°io) ® diag(mi, ■ ■ ■ ,m Nf ). By an appropriate counting the physical D5-brane charge of the 
orientifold which is equal to 1, one gets 

/ h= l -N k = [ h= l -N. k , f h=l-(2N + 2). 

By summing over the all a k contour, 

<f h = (2 E N k + 2N ) + 2 = (2N + 2)-N f , <f h = -(2N + 2), 
Jp \ k=1 J Jq 



h = -l 



±m.i 



10 



where h is an one form on the Riemann surface. 

The extremum condition of W (2.9) where the coefficient (2N — 2) is replaced with (2iV + 2) 
will give rise to 

rQ fP ' / pi r-mi\ 

(2N + 2) ( k -(2N-2N f + 2) & " E / +/ Cfe = 0. 

"'po "'Po i=1 Vpo Jpo / 

There exists a function on the Riemann surface with an (2iV + 2)-th order pole at Q, an 
(2N — 2Nf + 2)-th order zero at P and simple zeros at x — ±mj for each i = 1, 2, • • • , Nf. The 
function is simply related to the resolvent [37, 38] 

z{x) = ^Jx 4 W'(x) 2 - f(x) - x 2 W'(x). 

This function has an (2iV + 2)-th order pole at Q and at least zero-th order pole at P since f(x) 
is a polynomial of at most (2iV — 2)-th order. Therefore f(x) should contain a factor (x 2 — mf) 
for each i in order for z(x) to have a simple zero at x — ±rrii and should be at most 2Nf-th 
order in order for z(x) to satisfy an (2A^ — 2Nf + 2)-th order zero at P. Combining these two 
conditions and putting the proportional constant as A^f^ 2 27V/+4 , the spectral curve is 

\ i=l / i=l 

In the gauge theory side, for the Seiberg-Witten curve [47, 50, 57] for Sp(N) case at least one 
hypermultiplet of exactly zero mass (for example, rriNf = 0), the above matrix model result 
is exactly the same the one in [47, 50, 57] because the extra piece which is peculiar to Sp(N) 
case, the product of quark mass term vanishes. For nonzero quark mass, one can consider the 
following resolvent z(x) = \J (x 2 W'(x) + g) 2 — f(x) — (x 2 W'(x) + g) by including the constant 
term (x independent term) g inside the square. By identifying this g with A^_t, 2 Nf Ylf^rrii, 
one obtains the Seiberg-Witten curve. 
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